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20.  ABSTRACT  (Continued) 

-These  same  data  have  been  compared  with 
tion  data  from  the  Wideband  satellite, 
filiation  data  varies  with  perturbation 
as  does  pV-  With  realistic  propagation 
inferred  perturbation  levels  can  be  made 
However,  the  long-accepted  unity  separat 
spectral  indices  is  not  observed.  This 
shortcomings  in  the  theory  or  to  lack  of 
the  two  measurements . 


nearly  simultaneous  phase  scintilla- 
The  power-law  index  of  the  phase  scin- 
strength  in  exactly  the  same  manner 
model  parameters,  the  scintillation- 
to  match  those  measured  in-situ. 
ion  between  the  in-situ  and  phase 
discrepancy  is  attributed  either  to 
temporal /spatial  comparability  of 
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SUMMARY 


Radiowave  scintillation  has  been  a  continuing  concern  for  DNA  because 
of  its  adverse  effects  on  transionospheric  communications  in  both  natural 
and  nuclear-disturbed  environments.  The  Wideband  satellite  experiment 
was  designed  to,  and  already  has,  clarified  several  aspects  of  both  the 
morphology  and  the  propagation  channel  characteristics  of  ionospheric 
scintillation. 

One  of  the  specific  objectives  of  the  Wideband  experiment  has  been 
to  establish  the  connection  between  ionospheric  density  irregularities 
and  the  radiowave  scintillation  that  they  produce.  The  theory  predicts 
that  in  a  power- law  environment  the  phase  scintillation  spectrum  is  a 
direct  two-dimensional  mapping  of  the  three-dimensional  irregularity 
spectrum,  but  this  has  never  been  effectively  demonstrated.  With  this 
goal  in  mind,  an  analysis  of  in-situ  data  from  the  AE-E  satellite  has  been 
carried  out  in  a  manner  consistent  with  the  Wideband  analysis.  These  data 
provide  an  unambiguous  one- dimensional  measure  of  the  medium  scale  density 
irregularities  during  strongly  disturbed,  nighttime  equatorial  spread-F 
conditions . 

The  analysis  has  verified  the  power  law  form  of  the  in-situ  density 
spectrum,  with  an  average  spectral  index  close  to  that  measured  by  ether 
experimenters.  In  addition  to  this  expected  behavior,  however,  we  have 
found  an  apparently  new  result:  a  systematic  decrease  in  the  power- law 
slope  with  increasing  turbulence  level.  This  dependency  of  the  slope  on 
perturbation  level  of  the  in-situ  irregularity  spectrum  has  important 
ramifications  both  for  nuclear  phenomenology  and  propagation  prediction. 

A  similar  decline  of  phase  spectral  index  with  increasing  disturbance 
has  been  observed  for  some  time  in  the  Wideband  data,  but  had  previously 
been  attributed  to  diffraction  effects.  Now,  the  effect  seems  to  reflect 
a  true  in-situ  characteristic.  In  fact,  the  phase  spectral  indices  for 
the  same  period  of  time  as  the  AE-E  observations  show  almost  identical 
behavior  as  the  density  spectra. 
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By  using  a  spectral  approach  and  propagation  theory  based  on  a 
common  turbulent  strength  parameter,  the  one-dimensional  in-situ  and 
two-dimensional  phase  spectra  can  be  related.  In  practice,  however,  the 
very  different  orbits  of  the  two  satellites  make  such  comparisons  diffi¬ 
cult.  Nonetheless,  the  in-situ  turbulence  levels  are  consistent  with 
those  inferred  from  the  Wideband  data. 


PREFACE 


JSL 


The  authors  gratefully  acknowledge  the  assistance  of  Drs.  William 
Hanson  and  J.  P.  McClure  of  the  University  of  Texas  at  Dallas,  in  the 
selection  and  interpretation  of  the  AE-E  data. 

This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E 
RMSS  Code  B322080462  I25AAXHX64017  H2590D. 


3 


TABLE  OF  CONTENTS 


SUMMARY  . 

PREFACE  . 

LIST  OF  ILLUSTRATIONS  . 

I  INTRODUCTION  . 

II  IN- SITU  ANALYSIS  AND  RESULTS  .  .  . 
Ill  SCINTILLATION  ANALYSIS  AND  RESULTS 

IV  DISCUSSION  . 

REFERENCES  . 


4 


4 


LIST  OF  ILLUSTRATIONS 

1  Comparative  Spatial  Scales  of  AE-E  and  Wideband 

Measurements . . .  8 

2  Typical  AE-E  Pass,  Orbit  15058,  372-km  Altitude  .  10 

3  AE-E  Density  Spectra  for  Two  Different  Detrend 

Cutoff  Frequencies .  13 

4  Local  Time  Variation,  AE-E  Data .  14 

5  Distribution  of  Observed  Spectral  Slopes,  In-Situ 

Density  Data . 

6  Observed  Spectral  Slope  vs.  Spectral  Strength, 

In-Situ  Density  Data .  16 

7  Log  RMS  In-Situ  Density  Perturbation  vs.  In-Situ 

Spectral  Strength  .  17 

8  In-Situ  Measured  Three-Dimensional  Turbulent 

Strength .  20 

9  Wideband  Phase  Spectral  Slope  vs.  S ,,  Intensity 

Scintillation  Index  at  137  and  378  MHz .  23 

10  Distributions  of  Wideband  Phase  Scintillation  Data, 

137  MHz .  24 

11  Phase  Spectral  Slope  vs.  Spectral  Strength  T  at 

137  MHz .  25 

12  Scintillation  Measured  Three-Dimensional  Turbulent 

Strength .  27 

13  Comparison  of  AE-E  One-Dimensional  and  Wideband 

Two-Dimensional  Spectral  Slopes  vs.  Three-Dimensional 
Turbulent  Strength .  29 


5 


I  INTRODUCTION 


Electron  density  irregularities  in  the  ionosphere  are  observed  over 
a  vast  range  of  spatial  scales.  At  one  extreme  are  the  thousands-of- 
kilometer  scale  perturbations  associated  with  tidal  motions  and  planetary- 
scale  waves.  At  the  other  are  meter  scale  irregularities  presumably 
associated  with  diffusive  (ion  gyroradius)  or  electrostatic  (drift  wave) 
mechanisms.  Between  these  limits,  and  extending  over  at  least  three 
orders  of  magnitude  in  scale  size,  is  a  spatial  scale  regime  over  which 
the  irregularity  size  distribution  can  be  characterized  as  power  law 
[Rufenach,  1972;  Dyson  et  al.,  1974;  Phelps  and  Sagalyn,  1976;  Crane, 

1976;  McClure  et  al.,  1977;  Fremouw  et  al.,  1978]. 

In  this  report  we  consider  the  central  portion  of  this  irregularity 
power-law  regime--the  kilometer  through  hundreds-of -kilometer  portion, 
as  measured  by  two  quite  different  methods.  One  measurement  is  a  one¬ 
dimensional  sampling  of  the  irregularity  structure:  low  spatial  resolu¬ 
tion  in-situ  density  from  the  Atmospheric  Explorer-E  satellite  (AE-E). 

The  second  measurement  is  two-dimensional:  multifrequency  scintillation 
from  the  DNA  Wideband  satellite.  The  spectral  range  of  the  in-situ  data 
corresponds  to  scale  sizes  larger  than  a  few  kilometers,  which  overlaps 
the  low  frequency  end  of  the  Wideband  phase  scintillation  spectrum  as 
illustrated  in  Figure  1. 

Two  weeks  of  data  are  analyzed  from  a  period  for  which  AE-E  was  over 
Kwajalein  Atoll  in  the  equatorial  Pacific,  at  times  roughly  coincident 
with  the  Wideband  observations.  The  measurements  correspond  to  a  local 
time  and  season  when  the  irregularity  onset  and  evolution  are  dominated 
by  the  Rayleigh-Tay lor  instability  (sec  Ossakow  [1979]  and  the  references 
cited  therein). 

In  Section  II,  the  in-situ  data  are  first  characterized  in  terms  of 
the  spectral  density  function  (SDF)  of  density  fluctuation.  Using  a 
three-dimensional  spectral  model  and  generalized  geometry  and  anisotropy, 
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the  measured  perturbation  strengths  are  related  to  their  corresponding 
three-dimensional  turbulence  strengths.  In  Section  III,  the  same  ap¬ 
proach  is  taken  with  the  scintillation  data  so  that  the  results  of  two 
methods  of  measurement  can  be  compared. 

The  most  important  finding  from  these  data  is  a  systematic  variation 
of  the  power-law  spectral  index.  From  the  AE-E  data,  the  median  value 
of  the  spectral  index  is  slightly  less  than  2.  However,  the  spectral 
index  varies  monotonically  from  values  larger  than  2  to  values  signifi¬ 
cantly  smaller  than  2  as  the  perturbation  strength  increases.  A  similar 
trend  is  present  in  the  spectral  index  of  the  Wideband  satellite  phase 
scintillation  data. 


II  IN- SITU  ANALYSIS  AND  RESULTS 


During  late  July  and  early  August  1978  the  AE-E  satellite  was 
activated  over  Kvajalein  Island  (9°24’N,  167°28'E)  for  a  number  of 
nighttime  orbits.  The  data  that  follow  are  from  the  thirteen  of  those 
passes  that  encountered  significant  ionospheric  disturbance.  These 
observations  were  made  from  0800  UT  to  1400  UT  or  2000  LT  to  0200  LT . 

Among  the  instruments  aboard  AE-E  is  the  ion  drift  meter  described 
by  Hanson  and  Heelis  [1975],  from  which  angle-of-arrival  information  can 
be  gleaned  (for  example,  McClure  et  al.,  1977).  Also  available  from 
the  drift  meter  is  log  ion  dcnsity--the  data  we  have  used  in  this  study. 
The  spacecraft  system  loses  some  log  density  data  each  fraction  of  a 
second  during  drift  measurements.  This  loss  of  data  produces  a  compli¬ 
cated  pattern  of  missing  points;  however,  by  interpolating  over  a  span 
of  only  approximately  1/50  s,  it  is  possible  to  produce  a  record  of  log 
ion  density  with  an  effective  3-llz  sample  rate.  At  the  satellite  veloc¬ 
ity  of  ~  7.7  km/s,  this  velocity  corresponds  to  a  sample  every  ~  2.6  km 
along  the  satellite  path.  By  comparison,  the  high  spatial  resolution 
data  from  the  retarding  potential  analyzer  (e.g.,  McClure  et  al.,  1977) 
are  sampled  at  ~  30  m. 

Figure  2(a)  is  a  sample  of  the  log  density  data  over  Kwajaloin. 

For  this  particular  pass,  and  for  most  of  those  during  the  two-week 
period,  the  satellite  was  below  the  F-region  peak,  as  observed  on  iono- 
gram  records,  at  an  altitude  of  372  km.  This  record,  which  shows  the 
signatures  of  distinct,  isolated  plasma  bubbles,  is  typical  of  the  AE-E 
data  tor  the  more  strongly  disturbed  nights  during  the  Kwajaloin  scintil¬ 
lation  season.  Coordinated  measurements  between  AE-E  and  the  Kwajaloin 
radar  using  one  of  the  current  data  sets  has  shown  that  the  in- situ 
bubbles  correspond  to  I-m  backscatter  plumes  [Tsunoda,  1980]. 
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I  RAW  DATA  I 


FIGURE  2  TYPICAL  AE-E  PASS,  ORBIT  15058,  372-km  ALTITUDE,  (a)  Raw  log  density  and 
(b)  Depended  linear  density. 


Wo  wish  to  characterize  these  in-sit u  data  by  the  SDF  of  their 
density  perturbations.  Before  spectral  analysis  can  be  done,  however, 
it  is  necessary  to  separate  the  perturbation  components  from  the  slow 
changes  in  background  ion  density.  We  have  done  this  by  using  a  sharp- 
cutoff  filter,  a  method  which  lias  been  well  proven  in  the  analysis  of 
scintillation  data.  The  use  of  a  filter  has  the  advantage  over  polynomial 


ili-L  rending  in  Lhat  precisely  whore  in  frequency  and  how  much  the  SDF  is 


being  affected  by  L he  process  are  known.  In  practice  the  data  are  con¬ 
verted  to  Linear  electron  density  and  low-pass  filtered  with  subsequent 
subtraction  of  the  smoothed  coaiponent  from  t he  original  data. 

The  detrend  cutoff  frequency  must  be  low  enough  to  provide  an  ade¬ 
quate  range  for  determining  spectral  slope,  yet  high  enough  to  truly 
separate  density  trends  from  perturbations.  A  0.02-Hz  or  a  spatial  scale 
cutoff  of  roughly  400  km  has  been  used  in  the  routine  numerical  accumu¬ 
lations.  The  same  data  shown  in  Figure  2(a)  has  been  detrended  at  0.05 
Hz  in  Figure  2(b).  Subsequent  processing  consists  of  the  Fast  Fourier 
Transform  (FFT)  calculation  of  the  SDF  of  the  detrended  density  perturba¬ 
tions  and  extraction  of  the  spectral  strength  and  slope.  The  spectra  are 
calculated  over  a  time  span  of  about  85  s.  This  provides  adequate  spec¬ 
tral  resolution  between  spa  tia  1  scales  from  about  5  km  to  600  km,  and 
generally  assures  stationarity  at  the  detrend  filter  cutoffs  of  interest. 
In  accumulating  the  statistics,  each  span  overlapped  adjacent  spans. 
Finally,  spectral  fits  to  the  functional  form 

®(f)  =  T1f"pl  (1) 

are  made  by  calculating  the  least  squares  fit  to  the  smoothed,  log-log 
SDF  beLween  0.1  Hz  and  0.7  Hz,  which  is  well  above  the  detrend  cutoff, 
vet  below  the  noise  floor  frequency  for  a  disturbed  record.  Using  the 
least  squares  method  forces  consistency  in  a  large  number  of  spectral 
fits  that  cannot  be  obtained  by  eye. 

Of  concern  in  the  detrending  process  is  the  detrend  filter  response 
to  the  rapid  density  changes  between  the  bubble  walls  and  the  quiet  back¬ 
ground  layer.  Several  of  these  transition  regions  are  indicated  by  arrows 
in  Figure  2(b).  The  degree  of  ringing  at  these  transitions  depends  upon 
the  filler  cutoff  1 requency--becomi ng  worse  as  the  filter  narrows.  To 
verify  that  this  energy  does  not  contaminate  Lhe  in-situ  spectra  in  the 
■iope-lit  region,  a  lumber  of  spectra  at  various  filter  cutoffs  have 
be. -n  co  "piled.  An  example  of  those  comparisons,  for  a  data  span  in  which 


the  transitions  are  particularly  conspicuous,  is  shown  in  Figure  3  for 
cutoffs  of  0.01  Hz  and  0.05  Hz.  At  the  low-frequency  end  the  detrend 
cutoff  is  conspicuous,  but  otherwise  only  details  of  the  spectra  have 
changed.  In  the  fit  region  the  spectra  are  virtually  unchanged,  differ¬ 
ing  in  slope  fits  by  only  a  few  percent.  Thus,  we  feel  that  the  effect 
of  abrupt  transitions  on  the  summary  statistics  are  negligible. 

Note  that  the  longer  period  detrend  spectrum  shows  no  tendency  to 
flatten  at  its  lowest  frequencies.  This  indicates  that  if  a  systematic 
outer-scale  structure  size  is  in  the  density  data,  it  is  well  beyond 
400  km  (0.02  Hz). 

The  mix  of  quiet  and  disturbed  regions  in  the  density  data  produces 
a  broad  range  of  spectral  strength  T^,  varying  over  more  than  60  dB.  To 
avoid  the  more  weakly  disturbed  data  for  which  the  spectral  slopes,  p  , 
may  decrease  by  intersection  with  the  noise  floor,  we  have  rejected  data 
for  which  T^  <  195  dB. 

The  local  time  at  Kwajalein  for  these  AE-E  passes  ranges  between 
2000  LT  and  0200  LT.  When  plotted  as  a  function  ot  time  ^Figure  4(a)], 
the  T^  values  show  a  systematic  decline  through  and  beyond  local  midnight 
This  systematic  decline  is  the  general  trend  expected  in  the  nighttime 
equatorial  ionosphere  from  scintillation  and  rocket  measurements.  Note, 
however,  that  these  disturbance  l.  .els  in  the  :  led ium-to-largo  scale  re¬ 
gime  are  still  strong,  at  0200  LT  despite  tile  general  disappearance  of 
1-m  backseat  ter  two  to  three  hours  earlier.  The  systematic  decline  of 
spectral  strength  with  Local  Line  raises  the  question  whether  it  is  due 
simply  to  the  decline  in  the  background  density.  In  Figure  4(b)  the  mis 

>  f/  « 

normalized  perturbation,  (Z.N/N)~  “ ,  is  plotted  versus  time,  and  also 

shows  a  systematic  decrease  throughout  the  evening.  Therefore,  we  can 
attribute  the  observed  decline  in  spectral  slope  to  a  true  decline  in 
Ling  degree  of  turbulence  over  the  observing  window. 

(he  second  parameter  that  characterizes  the  irregularity  spectra 
is  tin-  slope  Pj  _ Eq .  (  L ) j  for  which  the  distribution  is  shown  in  Figure  5 
the  distribution  of  slopes  centers  about  p^  =  1.9  and  is  skewed  somewhat 
toward  oven  shallower  values.  This  mean  value  is  vorv  close  to  that 


FREOUENCY  Hr 

FIGURE  3  AE-E  DENSITY  SPECTRA  FOR  TWO  DIFFERENT 
DETREND  CUTOFF  FREQUENCIES 
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FIGURE  4 


LOCAL  TIME  VARIATION,  AE-E  DATA,  (a)  Spectral  strength 
T]  and  (b)  rnts  ^N-'N 
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SPECTRAL  STRENGTH,  T,  <IB  mks  units 

FIGURE  6  OBSERVED  SPECTRAL  SLOPE  vs,  SPECTRAL  STRENGTH,  IN-SITU 
DENSITY  DATA 


interval.  In  our  case,  a  sampling  rate  of  3  Hz  implies  a  functional 

dependence  indicated  by  the  clashed  line  in  Figure  7,  which  differs  from 

the  observed  dependence.  On  the  other  hand,  the  high  correlation  between 

2,  1/2 

(£N)“/  ~  and  1^  reassures  the  validity  of  the  characterization  of  the 

spectra  in  terms  of  T ^  and  p^. 

Up  to  this  point  we  have  dealt  with  the  two  measurable  parameters, 

T ^  and  p ^ ,  which  characterize  the  temporal  SDF  of  the  in-situ  density 
perturbation.  We  now  wish  to  consider  the  conversion  of  those  data  to 
a  spatial  parameter  that  can  be  generally  compared  to  other  measures  of 
turbulence,  and,  in  particular,  to  the  scintillation  data  of  the  next 
sec  t ion . 


6 


SPECTRAL  STREMGTH.  T, 


FIGURE  7  LOG  RMS  IN-SITU  DENSITY  PERTURBATION  vs.  IN-SITU  SPECTRAL 
STRENGTH.  The  dashed  line  indicates  predicted  correlation  for  constant 
spectral  slope. 


In  Lite  tumoral  case,  the  time-space  conversion  of  the  in-situ  data 
involves  both  the  direction  of  the  spacecraft  and  the  irregularity 
anisotropy.  This  can  he  shown  starting  from  an  appropriate  functional 
form  of  the  three-dimensional  fluctuation  spectrum,  e.g.,  the  general 
form  for  the  power- law  SDK  in  three  dimensions  of  Kino  and  Fremouw  "1977 
Because  v.v  are  investigating  a  spatial  regime  well  above  the  inner  scale 
v.e  have  chosen  the  simplification  of  that  model  used  by  Kino  1979  . 

Kino  1979  retains  generalized  anisotropy,  and  geometry;  he  de¬ 
scribes  the  three-dimensional  fluctuation  SDK  as 

.N(q)  =  aht  q"  '+l  )  t  -’) 

when  no  outer  scale  cutoff  is  assumed.  The  terms,  a  and  b,  describe  the 
anisotropy  of  the  densitv  irregularities  and  are  the  axial  ratios  alone 
and  icross  the  magnetic  f i e  Id  ,  respectively. 

C  -  8. 3/2  q2v'Jr(v  +  L  /  J )  /  ~  (v  -  1) 

s  o 
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i  ..  th.  •,  t  r.  ii;',  l  li  oi  the  three-dimensional  in-situ  turbulence  dependent 
.ip. mi  ih.  i)uLi-r  sc. ili'  ki.'i  number,  q  ,  and  the  slope  parameter,  v. 

As  liine  l't/u  imli  cat  i’ll ,  the  thri-i'  dimensional  autocorrelation 
lumlion  lh.it  lollows  iron  Eq .  (d)  is 


It.  +  1/J) 


K.-l(qoy)  * 


ici.  i  t  ran  i  or.it  i  on  .a  Eq,  t  *)  with  the  spatial  variable  replaced  by 
1  oiit.-ii:  o  product  vie  Ms  tin  oiu  -d  i  men  s  i  ona  1  temporal  SDK  as  an  in- 

...  pi  o  Id  e.i  .tire  i  t  : 


i  i  „  -  I'd) 


t  1  ■  I  q~  +  i/v  )") 

P  l  o  p  / 


-  1/d  .  (5) 


i  s  >:■'  -i  il  l.!  spied,  ,  is  dependent  upon  both  thi'  probe  velocity  and 
I  '  I  I  I  :  I  :  I  ■  I  r  i  t  v  nil  .otr.'i'V.  1'ro:  i\i  ilo  .  I  ‘1 7u  j  ,  the  for;:,  of  v  is 

p 


l  h,  p.  •  l.'iilv  and  tin-  matrix,  C,  i derived  iron  till’ 

:  ill  d  i:  1  i ;  its  oil  ■[  that  Kino  and  Eremouv  L 1 1 7  7  d  i  scuss. 

ii!  : :  1 1  (o  i.  that  lor  rod-like,  i ield-a  llgned  irregularities 

i  ill  iii  ",  i,  tii,  oi  : -d  ia.-ona  1  tii'is  oi  t  are  aero.  Furthermore, 

i  i  opi  i  ii:  si  ii  at  i,  (a  -  1)  the  throe  di  a-,  tonal  terns  are  unity, 


i  i  pi.  th.  pro!'  ,  |»  i  ,  d  . 


iii  iib  ...  II-  o 


!  in  tlio  outer  scale  cuLum,  Eq .  (5) 
to  Eq .  i  1),  used  to  ch.iractri'i/.i'  the 


(  I  )  -  1,1 


-(dv-l) 


1 


(8) 


C  7(v  -  1/2) 

X  =  - 

1  9  ?v_l 

4.rr(v  +  1/2)  v  (2jc/v  )“V 

Note  from  Eq.  (7)  that  the  slope  parameters  relate  as 


v 


pi  +  1 


(9) 


Eqs.  (7)  and  (8)  concisely  relate  the  temporal  SDF  measurables,  1' 
and  p^,  to  the  three-dimensional  spatial  turbulence  level  Cor  generalized 
anisotropy  and  geometry.  In  applying  the  measured  spectral  character¬ 
istics  to  calculate  the  strength  of  turbulence  via  Eqs.  (7)  and  (8),  we 
have  chosen  to  use  the  measured  values  of  spectral  strength  with  a 
functional  dependence  of  slope,  p^(log  T  ),  based  on  the  least  squares 
fit  illustrated  in  Figure  0.  The  resulting  distribution  of  C  ,  based  on 
the  in-situ  measurements,  is  shown  in  Figure  8(a).  The  distribution  is 
broadened  from  tiiat  of  T  ,  primarily  from  the  nonconstant  slope  which 
has  been  used,  and  secondarily  from  the  variation  in  satellite  direction, 

and,  therefore,  v  from  orbit  to  orbit. 

P 

Mom'  instructive  is  the  dependence  of  spectral  slope,  p^,  on  C  as 

shown  in  Figure  8(b).  The  decrease  of  spectral  slope  with  increasing 

turbulence  strength  is  again  clear;  there  is  an  overall  decrease  of 

about  unitv  slope  between  2.3  and  1.3  over  about  t>0  dB  of  increase  in  C  . 

s 

Because  the  direction  of  AE-E  is  primarily  magnetic  west-east  at  Kwaja- 
Lein,  Figure  8(b)  changes  Little  with  the  choice  of  anisotropy  model. 

From  isotropic  to  1000:1  or  oven  nonrod  models,  the  C  values  change  by 
only  a  few  percent,  and  the  dependence  of  slope  upon  C  indicated  bv  the 
Least  squares  fit  remains  virtually  unchanged. 

To  summarize  the  in-situ  data,  we  observe  a  decrease  in  the  degree 
of  turbulence  with  time  over  the  evening  hours  from  2000  LT  to  0200  LT. 
The  spectral  index  of  the  density  pertu rbaL ions  show  an  arithmetic  mean 
very  close  to  that  observed  by  others  .  When  we  plot  the  spectral  index 
a  function  of  spectral  strength  of  three-dimensional  turbulence,  however, 
we  observe  a  systematic  decrease  in  spectral  slope  with  increasing 


turbulence. 
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Ill  SCINTILLATION  ANALYSIS  AND  RESULTS 


During  the  same  period  that  the  AE-E  data,  just  described,  were 
collected,  the  DMA  Wideband  system  [ Fremouw  et  ai.,  1 9 7 7 J  was  operating 
at  Kwajalein.  The  scintillation  season  at  Kwajalein  is  centered  about 
local  summer,  and  shows  a  considerable  degree  of  variability  from  week 
to  week.  During  the  late-JuLy /ear ly-August  L978  period  of  interest,  the 
overall  scintillation  was  moderate  to  strong.  Considerably  more  severe 
scintillation  has  been  observed  at  Kwajalein  during  other  periods 
^Livingston,  L980j.  Of  the  satellite  data  collected  daily,  nine  night¬ 
time  passes  were  disturbed  and  form  the  data  base  for  the  following  dis¬ 
cussion.  These  were  collected  between  about  1130  UT  and  1300  Jf  (.1300  LT 
to  0100  LT),  thereby  overlapping  the  AE-E  observa tions . 

Fremouw  et  al.  ,1978j  review  the  analysis  applied  to  the  Wideband 
data.  in  brief,  the  phase  coherent  system  allows  measurement  of  signal 
intensity  and  phase  at  several  VHP  through  L-band  frequencies.  A  sharp 
cutoff  filter  at  10  s  is  used  to  separate  the  rapid  phase  components 
fro-  the  slowly  changing  background.  This  scintillation  component  is 
su bsequcntly  processed  in  30-s  blocks.  Included  in  this  processing  is 
the  computation  of  phase  spectral  density  functions  at  137  Miiz  and 
37S  o'l,:,  which  are  Least  squares  fit  over  a  selected  frequency  range  to 
th.'  form 

rtf)  =  T  r"p  ,  (10) 

analogous  to  the  AE-E  processing  „  Kq .  tOJ- 

Kino  _  197-1  j  has  demonstrated  and  discussed  the  consistency  of  the 
phase  spectra  L  parameters  between  fill',  l' 111',  and  L-band  frequencies.  In 

) 

particular,  he  shoes  Lli.it  p  is  independent  of  frequence,  and  that  T  \“ 
slum  the  small  correction  for  finite  reference  frequency  is  included. 


Both  of  those  items  indicate  that  diffraction  effects  are  probably  negli¬ 
gible  in  the  phase  spectra. 

In  Figure  9  we  demonstrate  that  same  consistency  for  the  much 
smaller  current  data  set,  where  p  is  plotted  versus  the  S,  intensity 

•4 

.scintil lation  index.  Least  square  fits  to  the  data  are  indicated  for 
137  MHz  and  378  MHz.  In  general,  the  agreement  between  the  slopes  at 
the  two  frequencies  as  a  function  of  S,  is  excellent.  Given  this  con- 

-4 

sistency,  and  to  simplify  discussion,  we  will  present  only  the  137 -MHz 
data  in  what  follows.  Similar  analysis  for  the  378-MHz  phase  scintilla¬ 
tion  lias  produced  quantitatively  equivalent  results. 

In  Figure  9,  for  the  low  S,,  weakly  disturbed  data,  the  phase  slopes 
are  flattened  by  intersection  with  the  spectral  noise  floor.  At  the  high 
S(,  stronger  scatter  end,  the  phase  slopes  also  decline;  this  effect  is 
more  conspicuous  in  these  data  than  in  the  larger  data  set  analyzed  by 
Kino  1979..  Kino  1979  attributed  this  effect  to  diffraction-produced, 
i ap  id ly  changing  steep  phase  fronts,  which  drive  the  SDF  dependence 
towaid  f  .  In  light  of  the  AE-E  Jata,  however,  this  trend  appears  to 
reilect  a  consistent  characteristic  of  the  in-situ  irregularities. 

As  was  the  case  for  the  in-situ  data,  the  patchy  nature  of  the 
equatorial  disturbance  produces  a  spread  distribution  of  the  observed 
scintil lation  strength.  To  eliminate  the  weak  scatter  disturbances, 
we  have  sorted  the  data  according  to  S,  index,  eliminating  points  with 
th  0.9  for  which  the  derived  phase  spectral  slopes  may  be  biased  by 
noise.  With  this  criterion  applied,  we  get  the  distributions  of  spectral 
strength  T  anti  slope  p  illustrated  in  Figure  10.  The  mean  value  of 
slope  here  is  ~  3.9,  only  slightly  smaller  tlan  the  value  that  Kino  [1979] 
t omul .  However,  it  is  considerably  lower  than  the  nominal  p  =  3  that 
Crane  _l97b_  reported.  In  Figure  11  the  observed  phase  scintillation 
slopes  are  plotted  versus  spectral  strength  to  produce  a  counterpart  to 
Figure  i)  of  the  in-situ  data.  Here,  as  in  the  other  data  we  find  a 
systematic  decrease  in  spectral  slope  with  increasing  turbulence  level. 

file  conversion  of  these  data  to  three  dimensional  turbulence  strength, 
C  ,  provides  a  means  of  direct  comparison  to  those  derived  in  situ. 
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DEBAND  PHASE  SPECTRAL  SLOPE  vs.  Sd  INTENSITY  SCINTILLATION  INDEX  AT  137  and  378  MH/ 


SPECTRAL  STRENGTH,  T  --  dB  mks  units 

FIGURE  11  PHASE  SPECTRAL  SLOPE  vs,  SPECTRAL 
STRENGTH  T  AT  137  MH/ 


In  the  phase  scintillation  case,  the  jn-situ  spectral  form,  F.q .  (2)  can 
ho  related  to  the  SDP  of  phase  using  the  phase-screen  model,  as  done  hv 
Rino  1 7 1>  j .  Analytically,  an  appropriate  form  Lor  the  three-dimensional 
autocorrelation  i unction  (ACF)  is  estahl  ished ,  from  which  the  two- 
dimensionaL  Lemporal  ACF  is  derived.  The  spa t ia l / tempera  1  conversion 
involves  the  relative  drift  of  the  medium  across  Lite  propagation  path , 
from  source  mo Li on  or  irregularity  drift.  The  resulting  temporal  ACF 
of  observed  phase  is: 


-O') 


T 


UD 


where  I  is  tile  temporal  outer  scale  defined 
1,  tlie  spectral  strength,  is  defined  as 


1  =  (v  ,-rq  /2  0,  and 

o  ell  o 


I 


•;•)  t;t: 


~(v) 


2v-l 


(2  )Jv+1:'(v  +  1/2) 


•f  f 


(12) 


In  Eq .  (12).  r  is  the  classical  electron  radius,  X  is  the  radio  wave- 
e 

length,  I,  is  the  equivalent  layer  thickness,  and  9  is  the  zenith  angle 


25 


o‘  propa  ■>,  at  ion  Lhronnh  tha  Lavar.  C  and  a  ,u'i'  toowl  r ica  L  UicLors 

all 

Chat  c.i'i  ahanpa  thronphonl  a  Wiilcli.i'iJ  satoLLito  pass.  C  rot  Lih  Is  lha 
dislnrhanaa  anhanaamant  prod  load  hv  aohafanl  propapa L  i  on  ,  a.p.,  .ilu'is; 
tho  lonp  axis  of  an  oMoiulod  L  rrapu  Lari  lv ,  whi  la  v  __  aivou’its  la.'  l  ho 
spat  ia  l-l.anpora  L  amrorsion  at  Lha  ionosphari  a  ponol  ra  L  ion  paint.  i>  - 
oansa  v  dapands  upon  lha  propagation  d  i  roo  l  ivm  with  rospott  to  lha 
principal  i rrapu  larilv  axis,  it  is  also  highly  dopotulaiu  upon  th  ■  alia  s,  n 
an  isotropy  nodal. 

Kawrilinp  Kq.  (.  L  L )  lot  fraqaanc  ias  wo  L  L  ha  Low  Llia  out.  r  sra  la  :  i  « 
a  form  lor  tho  pltaso  SI)/  similar  to  Kq .  tL)  ; 

-  ’\ 

A  f )  =  i  f  .  I  I  ' ) 

tdvrparls  xi  with  Kq ,  (,l»)  shows  that  in  t  ha  phaso  so  i  iu  i  1  la  L  i  an  aa-i  , 

p  -  A.  j 

In  Kq .  (  IJ) ,  thoiv:  •»  v,  wo  haao  a  -noons  to  roncorl  tiia  .a.i  ;u:v.i  1  | 

and  p  La  t  hrao-d  imoiis  ioua  l  luronLanL  stivn.'.th  C  ,  as  doao  with  tha  in-  ,! 

i! 

situ  Jala .  As  in  Lha  in-si  in  analysis,  tho  aioasu.Wol  raluos  of  1'  ha.o  j 

boon  usad  to  aalatilato  a  usinp  an  axpi'assion  for  slopa,  p  t  La  *  f),  ila- 

s 

riood  l  ro:n  tha  l  oas  t  -  sqna  ros  it  of  i'ipnra  11.  For  initial  compa  r  i  sons 
wo  horn  assusa.l  an  aquiaaLant  .'»:>«)— k-n  lliiok,  uni  fo"i:i  lavar  oantal’ad  at 
>'.•  *  k;n,  and  10.1:1  i  n-or.u  l.ari  tv  rod  modal.  Noto,  I  ros  Kq .  that  C 

is  Linaarly  ihpainlont  upon  tho  lavar  thioknoss  L.  Mo  vovor,  fo r  tha  rad¬ 
ii  ka  matlal,  tha  i  rtv.vu  lari  lv  anisotropy  sfonpLy  at  foots  tho  C  aalcuLa- 

s 

lion  throu.O)  lha  r  Lana  onlv  at  low  axial  ratios.  liovotul  a  '  L0  tha 

i’ !  I 

all  act  saturalas  rapidly  lioauuso  o!  tha  north-souLh,  noarly  f ia Id -a  Li pnad 
di  faction  ai  satalliLa  travel.  If  wo  raslriaL  oursohvs  then,  to  rad¬ 
ii  k'.t  inv.’.u  l  aril  Los  and  realistic  Lavar  thicknesses  and  heiphls,  lha 
valno  al  i:  aan  ho  shil  tad  only  soma  law  till. 

The  phaso  sc int i I  La t ion-imp  1 ied  distrihntion  of  C  is  shown  in 
I’ipnra  Ldta).  Lt  is  of  ronpliLy  tha  samo  niapnilude  as  Lha  in-situ  imp  Li  oil 
C  data,  hnL  of  cpiiLa  different  shapo.  1L  ran  ho  arpnotl  that  tha  dif- 
taraiiL  disLrlhuLian  shapos  aro  consistent  with  Lha  pvoha  and  propagation 
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geometries.  The  probe,  travelling  west  to  east,  encounters  multiple 
plumes  varying  in  strength  and  internal  structure.  Over  the  thirteen 
orbits  used  here,  each  covering  over  3000  km,  an  approximately  uniform 
sampling  is  expected.  Wideband,  on  the  other  hand,  in  its  rapid  (15  min) 
north- to-south  scan  tends  to  sample  fewer  plumes,  and  when  it  does, 
probes  the  long  and  more  consistently  structured  axis  of  the  disturbance. 
The  conversion  from  T  to  C  also  produces  Figure  12(b),  in  which  the 
measured  phase  SDF  slopes  are  plotted  versus  C  .  The  decrease  of  p  with 
increasing  C  is  clear,  as  it  was  with  the  in-situ  data. 
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TV  DISCUSSION 

The  in-situ  and  phase  scintillation  perturbation  strengths  have 
been  expressed  in  terms  of  Cg,  a  common  measure  from  which  the  geometry, 
scale  size  and  anisotropy  dependences  have  been  removed.  In  Figure  13 
the  least  squares  estimates  and  their  rins  errors  from  the  in-situ  and 
phase  data  of  Figures  8  and  12  are  superimposed.  The  slopes  of  the 
least-squares  fits  p^(log  Cg)  and  p(log  C  )  differ  only  by  a  few  percent. 
The  agreement  between  these  results  from  two  very  different  experiments 
is  important.  It  firmly  demonstrates  that  a  constant  rate  of  decrease 
in  the  in-situ  spectral  index  with  increasing  turbulence  strength,  is  a 
consistent  property  of  medium-scale  spread-F  irregularities. 


FIGURE  13  COMPARISON  OF  AE-T  ONE-DIMENSIONAL  AND  WIDEBAND 

TWO-DIMENSIONAL  SPECTRAL  SLOPES  vs.  THREE-DIMENSIONAL 
TURBULENT  STRENGTH 
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The  theoretical  considerations  of  Sections  II  and  III  indicate  that 
the  in-situ  and  scintillation  phase  slopes  should  differ  by  unity.  We 
observe  instead  a  consistent  slope  behavior  for  both  measurements,  sepa¬ 
rated  by  about  0.6.  There  are  a  number  of  reasons  that  may  contribute 
to  this  discrepancy.  An  overall  decrease  in  the  phase  spectral  index 
may  occur  because  of  spatial  smearing  along  the  propagation  path. 
Diffraction  effects  in  phase,  ignored  here,  may  be  present.  Alternatively, 
the  three-dimensional  spectral  model  and  the  model  parameters  we  have 
chosen  can  be  questioned.  A  basic  step  in  the  long  accepted  theory 
relating  the  spatial  autocorrelation  function  of  density  to  that  of  phase 
is  that  the  structure  decorrelation  scale  is  small  compared  to  the  con¬ 
tributing  layer  thickness.  If  this  criterion  is  not  satisfied,  the 
assumption  of  unity  separation  of  the  phase  and  p^  +  2  in-situ  slope  may 
be  incorrect. 

It  might  also  be  argued  that  the  in-situ  and  scintillation  data 
sets,  in  terms  of  the  details  of  these  slopes,  are  for  some  reason  not 
generally  comparable.  This  could  be  a  result  of  the  very  different 
orbital  directions  of  the  two  satellites.  More  likely  it  may  be  due  to 
the  spatial  separations  in  latitude,  longitude  and  altitude  between  the 
regions  probed  by  the  two  experiments.  For  example,  although  the  in-situ 
and  phase  Cg  distributions  roughly  overlay,  it  is  conceivable  that  the 
scintillation  data  is  dominated  by  more  strongly  perturbed  and  more 
shallowly  sloped  structure  at  altitudes  not  sampled  by  AE-E. 

Resolution  of  these  spectral  slope  questions  may  lie  in  the  large 
volumes  of  scintillation  and  in-situ  data  which  already  exist.  Further¬ 
more,  systematic  comparisons  of  these  larger  data  sets  may  help  explain 
the  decrease  of  spectral  slope  with  the  increasing  level  of  turbulence. 
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